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ABSTRACT 

We report the detection of a non-zero time delay between radio emission measured by the VLBA at 15.4 GHz 
and 7-ray radiation (7-ray leads radio) registered by the Large Area Telescope (LAT) on board the Fermi 
Gamma-Ray Space Telescope for a sample of 183 radio and 7-ray bright active galactic nuclei (AGNs). For 
the correlation analysis we used 0.1 - 100 GeV 7-ray photon fluxes, taken from monthly binned measurements 
from the first Fermi LAT catalog, and 15.4 GHz radio flux densities from the MOJAVE VLBA program. The 
correlation is most pronounced if the core flux density is used, strongly indicating that the 7-ray emission is 
generated within the compact region of the 15 GHz VLBA core. Determining the Pearson's r and Kendall's 
T correlation coefficients for different time lags, we find that for the majority of sources the radio/7-ray delay 
ranges from 1 to 8 months in the observer's frame and peaks at approximately 1.2 months in the source's frame. 
We interpret the primary source of the time delay to be synchrotron opacity in the nuclear region. 
Subject headings: galaxies: active — galaxies: jets — radio continuum: galaxies — gamma rays: galaxies 
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1. EVITRODUCTION 

Long-term systematic observations of active galactic nuclei 
(AGN) have produced detailed light curves that are a power- 
ful tool for investigating the nature of these highly energetic 
phenomena. AGNs show variability across the full electro- 
magnetic spectrum and multi-frequency monitoring programs 
serve to establish connections between flux variations at dif- 
ferent energy bands. 

The first comparison between long-term records of radio 
(10.7 GHz) and optica l fluxes for a sample of 24 AGNs 
jPomphrev et alj 119761) showed a correlation in only 13 
sources, with the optical events preceding radio by intervals 
of to 14 months. Subsequent analysis of light curves for 18 
AGNs taken in the op tical and radio (4.6-14.5 GHz) domains 
(IClements et alJ [1995 ) demonstrated similar statistics: nine 
sources exhibited positive radio-optical cor relations, with a 
time lag ranging from to 14 months. iTornikoski et alj 
(|T994) reported on a correlation in 10 out of 22 sources com- 
paring optical and radio (22-230 GHz) observations. Re- 
markably, in 6 sources the variability was simultaneous. A 
more recent statistical study of the time delay between indi- 
vidual milli meter- and centimeter-wave flare peaks made by 
iHovatta eTa l. (2008) for a sample of 55 sources showed a 
large scatter of time lags that ranged up to several hundreds 
of days between 4.8 and 230 GHz, an d tended to i n crease 
with decreasing radio frequency band. iJorstad et alJ (120101) 
investigated the flaring behavior of the quasar 3C 454.3 and 
showed that optical outbursts led the 23 GHz flares by 15 - 
50 days, confirming an earlier result bv lRaiteri et alj (120081) . 
who found millimeter flux changes lagging behind the op- 
tical on a time scale of about 60 days. A longer delay of 
about ten months between optical and 37 GHz radio flux vari - 
ations was reported in this source by IVolvach eTaP (l2008h . 
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iLahteenmaki & Valtaoial (l2003h studied the connection be- 
tween 7-ray emission detected by EGRET and phases of cor- 
responding radio flares at 22 and 37 GHz. They reported that 
the highest levels of 7-ray emission were detected 30-70 days 
after the onset of the high-frequency radio flare. 

The Fermi Gamma-Ray Space Telescope, successfully 
launched in June 2008, has opened a new era in 7-ray astron- 
omy. A number of AGN radio/7-ray connections have been 
established on the basis of the first three months of Fermi sci- 
ence operations and quasi-simultaneous VLBI observations, 
namely that the 7-ray photo n flux correlates wit h the par- 
sec scale radio flux density ( Kovalev et al.l l2009h . and that 
the jets of the LAT-detected blaza rs have higher-than- average 
apparent speeds (median of 15 c; iLister et aL 2009a), larger 
appar ent opening angles (median of 20°; 'Pushkarev et al] 
20091). and higher va riability Doppler factors (mean of 20; 
Savolainen et alJ I2OIO1) . In addition, AGN jets tend to be 



found in a more active radio state within several months of the 
LAT- detection of their strong 7-ray emission (iKovalev et alJ 
120091) . A significant correlation was detected even between 
non-simultaneous measurement s of 0. 1 - 100 GeV photon flux 
and 8 GHz VLBA flux density (lKovalevli200 9). Nonetheless, 
there are several key questions that remain to be addressed, 
which include: (i) the dominant production mechanism(s) of 
7-ray photons, (ii) the exact location of the high-energy emis- 
sion, and (iii) the nature of the AGN 7-ray duty cycle. In this 
Letter we investigate the radio/7-ray time delays for parsec- 
scale jets and place constraints on the localization of the re- 
gion where most of the 7-ray photons are produced. 

Throughout this Letter, we use the term "core" as the appar- 
ent origin of AGN jets that commonly appears as the bright- 
est feature in V LBI images of blazars (e.g.. .Lobanov.. 19981 ; 
lMarscheiil2008h . 

2. THE RADIO DATA AND SOURCE SAMPLE 

The MOJAVE program (iLister et al.l l2009b') is a long-term 
VLBA project to study the total intensity structure changes 
and polarization evolution of extragalactic relativistic radio 
jets in the northern sky. The observed sources include a sta- 
tistically complete, flux-density limited sample of 135 AGNs 
(MOJAVE-1). AH the MOJAVE-1 sources have J2000 dec- 
lination 6 > -20° and a 15 GHz VLBA correlated flux den- 
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sity Scon- > 1.5 Jy (2 Jy for 6 < 0°) at any epoch between 
1994.0 and 2004.0. The weaker radio blazars (5corr > 0.2 Jy) 
detected by Fermi extend the complete MOJAVE-1 sample 
to MOJAVE-2. The monitoring list currently consists of 293 
sources, 1 86 of which ar e members of the First Fermi LAT 
catalog (IFGL, Abdo et a l. 2010a) that are positionally asso- 
ciated with AGNs. We note that the 186 sources do not repre- 
sent a complete sample selected on either parsec-scale radio 
flux density or 7-ray photon flux. 

3. RESULTS 
3.1. Radio/'j- ray delay in VLBI cores 

Apart from the median 7-ray photon and energy fluxes, 
the IFGL catalog provides flux history data, in the form of 
monthly binned 0.1 - 100 GeV photon flux measurements 
during the first 1 1 months of the Fermi scientific operations, 
which started on 2008 August 4. The time sampling of our 
VLBA radio observations is source-dependent: objects with 
more rapid structural changes (i.e., faster apparent speeds) 
are observed more frequently. There are only five sources 
in our sample that are monitored with a cadence more fre- 
quent than once every two months. Starting in early 2009, 
fifty five bright 7-ray detections (> lOcr) from the Fermi LAT 
3-mont h list positio nally asso ciated with bright radio-loud 
blazars (lAbdo et al.i r2009c a; Kovalevll2009h have been incre- 
mentally added to the MOJAVE program. More than half of 
these new LAT-detected sources have fewer than three epochs 
of radio observations during the Fermi era, which precludes 
the correlation analysis of individual light curves. Therefore, 
our study is based upon a statistical approach. 

Overall, we obtained 564 VLBA images and corresponding 
model fits for 183 bright 7-ray sources (Table[T]i within a pe- 
riod from June 2008 through March 2010. The parsec-scale 
structure, typically represented by a one-sided core-jet mor- 
phology, was fitted with t he procedure modelfit in the Difmap 
package (IShepherdlll997[) using a limited numbe r of primar- 
ily circu lar Gaussian components, as described bv lLister et alj 
(I2009bl) . 

We tested for possible correlations between the 7-ray pho- 
ton fluxes and 15 GHz VLBA core flux densities using the fol- 
lowing procedure: (i) we selected all pairs of measurements 
where the difference in the radio and 7-ray epochs lay within 
a restricted time interval, for instance, [-0.5, +0.5] month, 
where the negative sign indicates that the radio measurement 
precedes the 7-ray one, (ii) if more than one pair of fluxes was 
available for a source, we selected the one with the epoch dif- 
ference closest to the mean of the time interval. The procedure 
was then repeated iteratively by shifting the time interval by 
0.5 month each time. We performed a quantative analysis that 
confirmed that our data do not provide any bias towards pos- 
itive pairs of radio/7-ray epoch difference. We used a cutoff 
of SNR > 3 for the 7-ray photon flux measurements to avoid 
a bias due to the lack of sources that are both weak in radio 
and 7-rays, and to exclude the influence of low-quality data 
points. For each data set we calculated the Pearson's r and 
non-parametric Kendall's r correlation coefficients, together 
with a corresponding probability of a chance correlation (Ta- 
ble |2]l. A non-zero radio/7-ray time lag, clearly seen as a 
bump in the correlation versus delay curves (Fig. [U left top 
panel), ranges from 1 to 8 months. The smooth fitted curves 
were obtained by applying a three point moving average. 

To test the robustness of this result, we estimated the uncer- 
tainty value of the correlation coefficients. Since both the ra- 



Table 1 

Sample of LAT-detected MOJAVE AGNs 



Source 


Alias 


Fermi name 


z 


Redshift reference 


(1) 


(2) 


(3) 


(4) 


(5) 



0003 + 380 
0015-054 
0048-071 
0048-097 
0059 + 581 
0106 + 013 
0106 + 678 

0109 + 224 

0110 + 318 
0113-118 



OB -082 



4C +01.02 
4C +67.04 

4C +31.03 



IFGL J0005 
1FGLJ0017 
IFGL J0051 
1FGLJ0050 
1FGLJ0102 
1FGLJ0108 
IFGLJOllO 
1FGLJ0112 
1FGLJ0112 
1FGLJ0115 



7 + 3815 
4-0510 
.1-0649 
6-0928 

8 + 5827 
.6 + 0135 
+ 6806 
.0 + 2247 

+ 3207 
.5-1132 



0.229 Abdo et al. (2010c} 
0.227 Abdo et al. (2010c) 
1.975 Abdo et al. (2010c') 



0.644 
2.099 
0.290 
0.265 
0.603 
0.672 



Abdo et al. (2010c) 
Abdo et al. (2010c^ 
Abdoetal. (2010c} 
Abdo et al. (2010c} 
Abdoetal. (2010c} 
Abdo et al. (2010c5 



Note. — Columns are as follows: (1) lAU name (1950); (2) alternate name; (3) IFGL 
name; (4) redshift; (5) literature reference for redshift. TablefTlis published in its entirety in 
the electronic version of the Astrophysical Journal. A portion is shown here for guidance. 



dio flux densities and the 7-ray photon fluxes are far from be- 
ing normally distributed, direct methods like the Fisher trans- 
formation or Student's t-distribution could not be used. Thus, 
we applied randomization techniques based on permutation 
tests to construct confidence intervals on the correlation coef- 
ficients. For each data set the randomization was done in the 
following manner: (i) we randomly swapped the radio mea- 
surements for one source with another source, keeping the 
7-ray fluxes the same; (ii) we calculated correlation coeffi- 
cients r and r from the randomized data. We then repeated 
these steps 2000 times. A 95% confidence interval (given in 
Table lU for the correlation coefficients was defined as the in- 
terval spanning from the 2.5-th to the 97.5-th percentile of the 
re-sampled r and r values. To estimate the null-basis level 
of the flux-flux correlation (whi ch is present due to an over- 
all radio/7-ray correlation; see § 13. 2t , we shuffled the 7-ray 
photon fluxes among the 1 1 measurements available for ev- 
ery source, keeping the epoch dates and radio flux densities 
the same. The resulting values were ro = 0.37 and tq = 0.26 
(Fig.E left bottom panel). When we randomly selected 90% 
(and 80%) of the sample the correlations remained signifi- 
cant, indicating that they are not driven by outliers. Addition- 
ally, we found no significant correlation between redshift and 
the VLBA core flux density averaged for the sample over the 
Fermi era. 

The wide range of delays in which the flux-flux correla- 
tions are significant (Fig [T] left) is presumably a result of the 
multiple parameters that determine the conditions in the nu- 
cleus (black hole mass, its spin and accretion rate), and in the 
nearby interstellar medium. The delay is also affected by ge- 
ometry, including the angle of the jets to our line of sight and 
the wide range of redshifts in our sample. The redshifts are 
known for more than 90% of the sources (166 out of 183). 
We re-did the analysis in the source's frame by dividing the 
radio/7-ray epoch time difference for each source by a factor 
of (1 +z). This gave a tw)ical time delay of ~ 1.2 months in 
the source frame (Fig. [TJ right), which corresponds to ^ 2.5 
months (for z ^ 1) in the observer's frame. The other sub- 
peaks are not significantly different from the null level of cor- 
relation, though they may indicate longer delays in a smaller 
number of sources. Note that the points on the radio/7-ray 
correlation curve are dependent. We expect the smearing even 
in the source frame delay because the core size is redshift 
dependent: more distant sources have higher rest-frame fre- 
quencies and, therefore should have a smaller core radius and 
shorter rest-frame delay. This effect can potentially be stud- 
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Figure 1. Left: Upper panel: flux-flux correlation level as a function of the time interval between integrated 0. 1 — 100 GeV photon flux and radio core flux density 
measurements in the observer's frame. Filled circles represent the Pearson's r statistic, and open circles denote the Kendall's r statistic. The error bars are shown 
at the 68% level, while shaded areas represent 95% confidence regions. The thick curves are constructed by applying a three point moving average. Dashed lines 
show the null-basis level of the coiTelation due to the overall 7-ray photon flux versus radio flux density correlation. Lower panel: estimation of the null level 
correlation obtained by shuffling 7-ray photon fluxes for every source and keeping radio flux densities the same. Right: same curves, corrected to the source rest 
frame. A sharp peak at a delay of ~ 1.2 months is detected. 



led in more detail using a larger sample subdivided into high- 
and low-redshift bins. Since a distribution of different delay 
values is expected, we do not estimate an error range on the 
detected peak. The peak value should only be taken as a typi- 
cal one for the AGNs in our sample. 

In Figure |2] we plot the integrated 0.1 - 100 GeV 7-ray 
photon flux against the 15 GHz VLBA core flux density 
for data pairs in which the 7-ray measurement leads the 
radio measurement by 2.5 ± 0.2 months in the observer's 
frame. The formal probability of a chance correlation is 
5 X 10"^. If we drop the points with photon flux greater than 
4 X 10"^ ph cm~^ s~', the correlation is still present at a very 
high level of significance. We also note the span over one or- 
der of magnitude in the 7-ray fluxes for a given radio flux den- 
sity. Several factors can contribute to this scatter: the shape of 
the 7-ray spectrum, K-correction effects, different seed pho- 
tons, and diffe rent Doppler boosting levels in the 7-ray and 
radio domains (lListerll2()07h . 

3.2. Localization of the ^-ray emission region 

Although radio loud AGN are generally known to have twin 
jet structures, the sources in our sample typically hav e a one- 
sided parsec-scale morphology jLister et alj l2009bl) . imply- 
ing strong selection effects and Doppler boosting of the jet 
emission. Localization of the 7-ray emission region within 
the AGN radio structure and its physical production mecha- 
nism(s) remain topics of active debate. There are three pos- 
sibilities for the site of the high-energy emission: (a) only 
within the unresolved radio core, (b) only in the resolved 
(downstream) jet region, (c) both in th e co re and jet. 

Repeating the same analysis as in § l3.1l for the total VLBA 
flux density, we found that the corresponding correlation co- 



efficients agree within the 95% errors with those found for 
the core flux densities. Thus, no firm conclusions can be 
drawn out of this comparison. Indeed, the sources in our 
sample are highly core dominated, with a median value of 
'Score/'SvLBA = 0-71, whcrc 5core and 5vLBA are the core and to- 
tal VLBA flux density, respectively. By contrast, when the jet 
flux densities (5vlba- 'Score) are used, the correlation coeffi- 
cients are significantly lower, making scenario (b) less prob- 
able. To reduce the uncertainty in the jet flux density esti- 
mations we excluded sources with a high (> 0.9) core domi- 
nance. Nevertheless, the correlation between the jet flux den- 
sities measured quasi-simultaneously with the 7-ray photon 
flux is still significant. This could be the result of both ra- 
dio and 7-ray regimes being boosted by similar beaming fac- 
tors (Lister et al. 2009a; Kovalev et al. 2009; Savolainen et al] 
I2OIOI) . Assuming the same Doppler factor for the jet and the 
core (opaque jet base), we would expect a weak but signif- 
icant correlation between the jet flux density and integral 7- 
ray emission. Alternatively, it is possible that for some low- 
redshift sources, 7-ray emission may also be pro duced out- 
side th e radio core, as could be the case in 3C84(Abdo et alJ 
l20091j) . where radio flare accompanying the 7-ray activity 
was detected in t he innerrn ost jet region. 

We note that iKovalevI (12009) found a correlation be- 
tween the Fermi LAT 0.1-100 GeV photon flux from the 
three-month integration and 8 GHz radio flux density non- 
simultaneously measured by the VLBA. The presence of this 
correlation in this highly variable population further suggests 
that Doppler beaming is the likely cause, and that the Doppler 
factors of the individual jets are not changing substantially 
over time. 
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Table 2 

Radio/7-ray flux-flux correlation statistics for tlie VLBA 
cores in tlie source frame 





N 


r 


T 


P 


(1) 


(2) 


(3) 


(4) 


(5) 



-4.9 ±0.2 
-4.5 ±0.2 
-4.0 ±0.2 
-3.5±0.1 
-3.0 ±0.2 
-2.5 ±0.1 
-2.0 ±0.2 
-1.5±0.1 
-1.0±0.2 
-0.5 ±0.1 
0.0 ±0.2 
0.5±0.1 
1.0±0.1 
1.2±0.1 
1.5±0.1 
2.0±0.1 
2.5±0.1 
3.0±0.1 
3.5±0.1 
4.0 ±0.2 
4.5±0.1 
5.5±0.1 
6.0 ±0.2 
6.5 ±0.2 
7.0 ±0.2 
7.5 ±0.2 
8.0 ±0.2 
8.5 ±0.2 
9.0 ±0.2 
10.0 ±0.2 
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27+0.02 
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24+^-^^ 
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1.6 X 10-' 

1.4 X 10"- 
2.0 X 10"' 

4.5 X 10-3 
1.9 X 10-2 
9.4 X 10-3 

3.6 X 10-2 

8.0 X 10-3 

1.1 X 10-2 
7.9 X 10-3 
9.6 X 10-* 



3.5 
7.3 
4.0 
1.6 
9.0 
1. 

2.6 



X 10-* 
X 10-^ 
X 10-^ 
X 10-*^ 
X 10"* 
X 10"* 
X 10-^ 



1.3 X 10-3 

2.4 X 10-3 

4.1 X 10-5 

4.4 X 10-5 
2.6 X 10-* 

2.2 X 10-3 

5.8 X 10-* 
1.6 X 10-3 

6.0 X 10-3 

1.5 X 10-2 

8.9 X 10-2 

2.1 X 10-2 



Note. — Columns are as follows: (1) radio-gamma time delay 
bin in months; (2) number of AGNs included in the bin; (3) Pear- 
son's r with uncertainties at a 95% confidence level; (4) Kendall's 
r with uncertainties at a 95% confidence level; (5) probabifity of 
a chance conelation. 



4. DISCUSSION 

What is the main source of the detected time lag? Disper- 
sion in the intervening medium cannot be responsible for a 
delay of several months. Our calculations following the Crab 
pulsar measurement by Abdo et al. (2010b) show that for an 
object at a redshift z = 1 the delay should not exceed a few 
seconds at 15 GHz. 

The most likely source of the observed time lag is syn- 
chrotron opacity in the nuclear region. The radio core is 
optically thick to synchrotr on emission up to the freq uency- 
dependent radius oc z/" (iBlandford & K6nigll [T979). This 
means that the 7-ray production region/zone is located up- 
stream with respect to the 15 GHz apparent core position. Be- 
ing induced by the same disturbance at a distance r-y from the 
black hole, the radio and 7-ray emission at their peaks are not 
observed simultaneously due to the opacity effect. Although 
the 7-ray photons escape immediately, it takes several more 
months for the perturbation to propagate farther along the jet 
until it reaches the r ~ 1 surface at 15 GHz radio emission 
(the radio core), and becomes detectable at radio frequencies. 
The corresponding distance travelled along the jet between 
the place where a 7-ray photon was emitted (r^) and the ra- 
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Figure 2. Integrated 0.1-100 GeV Fermi LAT photon flux versus 15 GHz 
VLBA core flux density for data pairs in which the VLBA flux density mea- 
surement was taken 2.5 ± 0.2 months after the LAT flux measurements. One 
point is plotted per AGN. 



dius of the radio core (rc) is 



ST/3cAt^_ 
(l+z) 



(1) 



where 6 is the Doppler factor, F is the Lorentz factor, /3c is 
the intrinsic jet speed, and Af*^^ is the observed time delay. 
Using the r elations for Dopple r factor S and apparent angular 
speed /3app dCohen et al.ll2007b , the expression (1) can also be 
written in the form 



Ar = 



„ A fsour 
app<-^'R-7 



smt 



(2) 



where /J^pp is the apparent jet speed, Af^°" is the radio to 
7-ray time delay in the source's frame, and 6 is the view- 
ing angle. Since we cannot measure radio/7-ray time de- 
lays for the sources individually from our data, let us con- 
sider a source with a set of parameters typical for a LAT- 
detected radio-loud blazar: appai-ent jet s peed /3app ^15 
(Lister et al."2009a"), viewing angle 6 ^ 3.6° (Push karev et al J 
2009; Hovattaetal. 2009), and a typical time lag of 1.2 
months in the source's frame. Under these assumptions we 
obtain the distance between the 7-ray production zone and 
the r « 1 surface at 15 GHz, the radio core, Ar ^ 7 pc, which 
corresponds to a projected distance of ^ 0.9 pc, or ^ 0. 1 mas 
for a source at z ^ 1. These estimates are consistent with 
the core radius obtained from the frequency dependent core 
shift measurements (Loba nov .1998; .Kovalev et al. 2008a,bl: 
lO'SulUvan & Gabuzdafl2009l) . This indicates that VLB! ob- 
servations at higher frequencies, for instance at 43 GHz 
and 86 GHz, should register shorter delay or even quasi- 
simultaneous flux variations with 7-ray flux at least for some 
sources, and might even resolve the region of the jet where 
7-ray emission is generated. 

The fact that the radio emitting region may be more phys- 
ically extended than the 7-ray emission zone may also con- 
tribute to the detected delay. Finally, we note that our analysis 
is sensitive to the delay between the peaks at radio and 7- 
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ray frequencies, but cannot address which flare, low- or high- 
energy, originates first. The latter analysis requires very well- 
sampled light curves for a large set of 7-ray-detected AGN. 

5. SUMMARY 

We have investigated the dependence between the inte- 
grated 0.1-100 GeV 7-ray photon flux and 15 GHz radio 
flux densities for a large sample of 183 LAT-detected AGNs 
observed by the VLBA within the MOJAVE program and con- 
clude the following. 

1 . The correlation between 7-ray photon flux and radio flux 
density is found to be highly significant. The correlation anal- 
ysis results for the core and total VLBA flux density are in- 
distinguishable. The correlation is systematically weaker if 
the jet flux density is used, providing further support for the 
localization of the 7-ray emission to the core region as well 
as a connection between the inverse-Compton 7-rays and the 
synchrotron radio emission from the jet. 

2. We found a non-zero radio/7-ray delay (7-ray leads the 
radio emission) that ranges from 1 to 8 months in the ob- 
server's frame and peaks at ~ 1 .2 months in the source frame. 
The delay is most likely connected with synchrotron opac- 
ity in the core region, although other mechanisms may play a 
role. 

3. The region where most of 7-ray photons are produced is 
found to be located within the compact opaque parsec-scale 
core. The typical distance between the 7-ray production re- 
gion and the 15 GHz radio core is estimated to be ^ 7 pc, 
which is consistent with the typical core radius derived from 
frequency dependent core shift measurements. 

These results a r e cons isten t with earlier fin dings reported 
by iKovalev et all (l2009h and iTavecchio et"an dlOlQi) which 
place further constraints on the localization of the 7-ray pro- 
duction zone in parsec-scale jets. The MOJAVE program is 
continuing to monitor a large sample of LAT-detected AGNs 
for use in more comprehensive studies associated with the 
next Fermi data release. 
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